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rj'(rj) — > 77 : A Tale of Two Anomalies 
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The radiative decays of the pseudoscalar mesons P = (7r° ,77,77') are of special interest as they provide an 
experimental window on the electromagnetic and colour chiral Ua(1) anomalies. tt° — > 77 is well described by 
the electromagnetic J7a(1) anomaly under the assumption that 7r° is a Goldstone boson for spontaneously broken 
chiral symmetry, but the analogous results for 7/(17) — * 77 are complicated by the colour Ua(1) anomaly in QCD. 
This paper reviews the theory of 77' (77) — > 77 decays, emphasising the role of the colour Ua(I-) anomaly and the 
renormalisation group. The relation to the Witten-Veneziano mass formula for the 77' and the QCD topological 
susceptibility is derived. The implications for the phenomenological analysis of P — » 77 are described and a 
proposal is made for a comprehensive re-analysis of existing data on P — » 77, P — > V7, where V = (p,oj,(j)) is 
a light 1~ meson, 77' (77) — > -KTvy, tjj — -> 77' (77)7, etc. Other applications, to the muon anomalous magnetic moment 
<7 M — 2 and the determination of the photon structure function gj in polarised deep inelastic scattering, are also 
briefly discussed. 
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1. Principal Results 

The radiative decays of the pseudoscalar 
mesons P — > 77, where P — (7r°, 77, 77') have long 
been of special interest because of their intimate 
relation with the chiral (axial) Ua(1) anomaly. 
The decay ir° — > 77 is of course the textbook 
example of the phenomenological importance of 
the electromagnetic contribution to the Ua(1) 
anomaly. It has played an important role both 
in understanding the theory of anomalies and in 
pinning down the quantum numbers of the quarks 
and number of colours N c in the early days of 
QCD. The decay rj' — > 77 is doubly interesting, 
since it also involves in an essential way the colour 
contribution to the Ua{^) anomaly. 

An immediate consequence of the colour J7a(1) 
anomaly in QCD is that even in the limit of mass- 
less quarks, the 77, unlike the 7r° and 77, is not a 
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Goldstone boson of spontaneously broken chiral 
symmetry. The success of the theoretical predic- 
tion relating the decay amplitude for ir° — ► 77 to 
the electromagnetic ?7a(1) anomaly is based on 
the pion being an approximate Goldstone boson 
in the massless limit. The complete theoretical 
analysis can then be abbreviated in 'soundbite' 
form by writing the familiar PCAC relation: 



(1) 



where onthe r.h.s. 7r is a phenomenological field 
for the 7T°. The corresponding analysis for the 
7/ is complicated first by the importance of non- 
trivial 77 — 77' mixing due to the explicit SU(3) 
breaking induced by m s m u ,md and second, 
by the presence of the colour contribution to 
the J7a(1) anomaly. Nevertheless, the analogous 
analysis can be pushed through using the anoma- 
lous Ward identities and, subject to reasonable 
dynamical assumptions, results in the following 
generalised PCAC relation for the 77 and rj': 



f^rnlptf +2n f G5 a0 



(2) 
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where 7j a = (77,77') and a = 0, 8. The new term 
on the r.h.s. reflects the presence of the anomaly. 
Very roughly, what happens is that due to the 
anomaly, the flavour singlet pseudoscalar mixes 
with the gluonic operator G^G^. The conven- 
tional PCAC relation applies only to the fictitious 
state before mixing. Rearranging this into a part 
involving the physical 77' therefore leaves a resid- 
ual 'glue' contribution, represented by G. How- 
ever, great care must be taken with this picture 
and especially with the interpretation of eq.(2). 
In particular, G is not to be understood neces- 
sarily as the lightest 0~ glueball, and the rela- 
tion (2) is only valid as an operator relation for 
insertions into zero-momentum Green functions. 
These points will be explained carefully in due 
course. 

In this paper we review the special features 
that arise due to the gluonic contribution to the 
Ua(1) anomaly when PCAC methods (which in- 
clude chiral Lagrangians) are used to analyse ra- 
diative pseudoscalar meson decays in the flavour 
singlet channel. In |l| , we presented an analysis of 
rf — > -y-y decay in the chiral limit of QCD, taking 
into account the gluonic anomaly and the associ- 
ated anomalous scaling implied by the renormal- 
isation group. Here, we summarise the results of 
a more recent analysis^] extending this to QCD 
with massive quarks, incorporating 77 — 77' mix- 
ing. In particular, we show how a combination 
of the radiative decay formula and a generalisa- 
tion of the Witten-Veneziano mass formula for 
the 77' could be used to measure the gluon topo- 
logical susceptibility x(0) in full QCD with mas- 
sive quarks. 

Our main result is summarised in the formulae: 

o. 

f aC " S»;°77 + 2n f A ffG 7 7 <*a0 = < m - (3) 

7T 

which describes the radiative decays, and 
f aa {m 2 ) aP f T ^=(2n f fA S a0 S b0 

(m u (uu) \ 

m d {dd) 
m s (ss) J 

(4) 

which defines the decay constants appearing in 
(3) through a modification of the Dashen, Gell- 



Mann Oakes Renner 2 ||,|| formula to include the 
gluon contribution to the Ua(1) anomaly. 

In these formulae, rj a denotes the neutral pseu- 
doscalars 7r°, 77, 77'. The (diagonal) mass matrix 
is (m 2 ) a /3 and j,o 77 is the appropriate coupling, 
defined as usual from the decay amplitude by 

(HIV") = -i9v a ll e *pa{lPlP2 eX (Pl) eP (P2)- The 

constant a^ m is the coefficient of the electromag- 
netic contribution to the axial current anomaly: 

d"J? 5 = d acb m c 4 + 2n f S aQ Q + a a cm ^F^F^(5) 

Here, J^ 5 is the axial current, 0!? = (775 T a q is 
the quark pseudoscalar and Q = ^trG^G^ is 
the gluon topological charge. m a are the quark 
masses (see eq.(16)). a — 0,3, 8 is the flavour 
index, T 3 ' 8 are ST/ (3) generators and T° = 1. 
The d-symbols are defined by {T a , T b } — d abc T c . 
Since this includes the flavour singlet Ua(1) gen- 
erator, they are only symmetric on the first two 
indices. For nf — 3, the explicit values are 

^000 = ^033 = ^088 = 2, C?330 = ds80 = 1/3, ^338 = 
^383 = —^888 = 1/V3. 

These formulae become clearer if we immedi- 
ately make the approximation m Ul md -C m s . 
The formulae for the n° decouple: 

7T 

where a 3 m = ^N c and 

f£ml = -(m u {uu) + m d (dd)) (7) 
leaving the following new results Q for the 77 and 

(8) 

where a° m = |iV c and a 8 m = -^=N C , together 
with the generalised Dashen formulae: 

(/°"') 2 "4 + (f° V ) 2m n = -4m s (ss) + 36A 

2 Hereafter simply referred to as the Dashen formula. 
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(9) 



The decay constants f aa in (3) are defined so 
as to satisfy the relation (4). It is crucial to 
recognise that in general they are not the cou- 
plings of the pseudoscalar mesons to the axial 
current Q. In the flavour singlet sector, such a 
definition would give a RG non-invariant decay 
constant which would not coincide with those in 
the correct decay formula (3). In contrast, all 
the quantities in the formulae (3), (4) are sepa- 
rately RG invariant g,g. The proof is not im- 
mediately obvious, and depends on the RGEs for 
the various Green functions and vertices defining 
the terms in (3), (4) being evaluated on-shell or at 
zero-momentum. 

As we have seen, since flavour SU(2) symmetry 
is almost exact, the relations for ir° decouple and 
are simply the standard ones with / 37r identified 
as f v . In the octet-singlet sector, however, there 
is mixing and the decay constants form a 2 x 2 
matrix: 



/ fOi]' fOr] 

faa I J J 

J -\ f 8 V > f8 r, 



(10) 



The four components are independent. In par- 
ticular, for broken SU(3), there is no reason to 
express f aa as a diagonal matrix times an orthog- 
onal rj — rf mixing matrix, which would give just 
three parameters. Several convenient parametri- 
sations may be made, e.g. involving two constants 
and two mixing angles, but this does not reflect 
any special dynamics. We return to this point 
when we discuss phenomenology in section 3. 

The novelty of our results of course lies in the 
extra terms arising in (3) and (4) due to the glu- 
onic contribution to the Ua(X) anomaly. The co- 
efficient A is the non-perturbative number which 
specifies the topological susceptibility in full QCD 
with massive dynamical quarks. The topological 
susceptibility is defined as 



X(0)= / d 4 x i(Q\T Q(x) Q(0)|0) 



(11) 



The anomalous crural Ward identities determine 
its dependence on the quark masses and conden- 



sates up to an undetermined parameter, viz. 



Aqq}J 

Notice how this satisfies the well-known result 
that x(0) vanishes if any quark mass is set to 
zero. 

The modified flavour singlet Dashen formula is 
in fact a generalisation of the Witten-Veneziano 
mass formula for the if. Here, however, we do 
not impose the leading order in l/N c approxi- 
mation that produces the Witten-Veneziano for- 
mula. Recall that this states 



2m K ^ ~ix(0)\ YM 

J 7T 



(13) 



To recover (13) from our result (see the first 
of eqs.(9)) the condensate m s (ss) is replaced by 
the term proportional to f^m 2 K using a standard 
Dashen equation, and the singlet decay constants 
are set to y / 2n// 7r . The identification of the large 
N c limit of the coefficient A with the non-zero 
topological susceptibility of pure Yang-Mills the- 
ory follows from large N c counting rules and is 
explained in ref.Q. Lattice evaluations || (see 
also ref.jTj]) in the quenched approximation find 
the following value, 



X(0) 



YM 



-(180MeV) 4 



(14) 



A similar result has also been obtained using 
QCD spectral sum rules || . 

The final element in (3) is the extra 'coupling' 
gc-y-y in the flavour singlet decay formula, which 
arises because even in the chiral limit the n' is 
not a Goldstone boson because of the gluonic 
[/4(f) anomaly. A priori, this is not the cou- 
pling of a physical particle, although (suitably 
normalised) it could be modelled as the coupling 
of the lightest predominantly glueball state mix- 
ing with if . However, this interpretation would 
probably stretch the basic dynamical assumptions 
underlying (3) too far, and is not necessary either 
in deriving or interpreting the formula. In fact, 
the gcyy term arises simply because in addition 
to the electromagnetic anomaly, the divergence of 
the axial current contains both the quark pseu- 
doscalar </>5 and the gluonic anomaly Q. Diago- 
nalising the propagator matrix for these operators 
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isolates the r\ and v( poles, whose couplings to 77 
give the usual terms g IJ77 and g v >yy However, the 
remaining operator (which we call G) also cou- 
ples to 77 and therefore also contributes to the 
decay formula, whether or not we assume that its 
propagator is dominated by a 'glueball' pole. We 
emphasise again that there is no need whatsoever 
to make any assumption about the spectrum of 
the G propagator in deriving the decay formulae. 

The presence of the coupling gcyy in (8) how- 
ever removes any immediate predictivity from the 
r/' — > 77 decay formula. In order to push the the- 
ory further, we therefore need to make additional 
dynamical assumptions. The obvious possibility 
is to explore whether the phenomenologically suc- 
cessful OZI rule can be applied in this context. 
(The OZI rule is of course closely related (but 
not identical^]) to the l/N c expansion, which is 
discussed in this context in ref. (!]].) In fact, as 
we now argue, the OZI rule does provide some 
justification for believing that the extra coupling 
5G77 is indeed small. The theoretical argument 
is based on the fact that gcyy is both OZI sup- 
pressed and renormalisation group (RG) invari- 
ant Since violations of the OZI rule are as- 
sociated with the Ua(^) anomaly, it is a plausi- 
ble conjecture that we can identify OZI-violating 
quantities by their dependence on the anomalous 
dimension associated with the non-trivial renor- 
malisation of J° 5 due to the anomaly. In this way, 
RG non-invariance can be used as a flag to indi- 
cate those quantities expected to show large OZI 
violations. If this conjecture is correct, then we 
would expect the OZI rule to be reasonably good 
for the RG invariant gcyy, which would there- 
fore be suppressed relative to <7»' 77 . (An impor- 
tant exception to this argument is of course the rj 
mass itself, which although obviously RG invari- 
ant is not zero in the chiral limit as it would be in 
the OZI limit of QCD.) Notice that this conjec- 
ture has been applied already with some success 
to the 'proton spi n' p roblem in polarised deep in- 
elastic scattering (lO). 

If it is indeed a good dynamical approxima- 
tion to assume gcyy is small compared to g n 'yy, 
we could combine eqs (8) and (9) to give a mea- 
surement of the non-perturbative coefficient A 
in x(0)- To see this, assume that the physical 



quantities , to,,/ , 3 ??77 and <jy 77 are all known 
and (temporarily) neglect <?g 77 . Clearly, the two 
purely octet formulae can be used to find f 8fl and 
f 71 if both g m -y and g^yy are known. The off- 
diagonal Dashen formula then expresses f 0ri in 
terms of f 0r ' . This leaves the two purely singlet 
formulae involving the still-undetermined decay 
constant f° n , the topological susceptibility co- 
efficient A, and the coupling gGyy The result 
follows immediately. If we neglect gcyy, we can 
find f° n from the singlet decay formula and thus 
determine A from the remaining flavour singlet 
Dashen formula. This is the generalisation of 
the Witten-Veneziano formula. Determining A 
in this way would of course be an important link 
between the phenomenology of 77' decays and the 
important subject of gluon topology in QCD ||. 

However, the issue of the magnitude of gcyy is 
ultimately an experimental question. It is there- 
fore crucial that phenomenological analyses of the 
data on r\ and rf decays start from the complete 
formulae (8), (9) and do a best fit involving the 
full set of parameters. Only then will we really 
know whether or not the extra anomaly-induced 
coupling gcyy is small. 

We return to this issue in section 3. First, we 
provide a PCAC-based proof of the key formulae 
(3), (4). Carefully used, these methods permit a 
very quick derivation. However, there are many 
subtleties in the analysis which are easily missed 
using these conventional techniques. For this rea- 
son, we prefer the more field-theoretic approach 
explained in refs.[QJ|] (see also ||). This makes 
very clear exactly what dynamical assumptions 
must be made and where they enter the argu- 
ment. The price is working with a formalism 
which, although elegant and easy to use, is un- 
familiar to many phenomenologists. We present 
this field-theoretic derivation in section 5, al- 
though readers interested only in the phenomeno- 
logical results will find the discussion in sections 
1-4 complete and self-contained. 
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2. U A (X) PCAC 

Consider first QCD by itself without the cou- 
pling to electromagnetism. The axial anomaly is 

d" J^ 5 = M ab 4 + 2n f Q5 a0 (15) 

The notation is defined in ref. Q . The quark mass 
matrix is written as m a T a , so 



= m°l + m 3 T 3 + m 8 T 8 (16) 



The condensates are written as 

(uu) 
(dd) 




This is achieved by 

G = Q-{Q<f%)({<hfa)& 6 (23) 
which reduces at zero momentum to 

G = Q + 2n f A$£4 (24) 



m u 











m d 











m s 



and we define 



_a rTaajf. — l ±b 



(25) 



With this choice, the (G G) propagator at zero 
momentum is 




where 
M ab = 



is the VEV (qT c q). Then 



dacbrrf 



(18) 



An essential approximation in the application 
of PCAC methods is to use identities valid at 
zero momentum and make certain smoothness 
assumptions (see section 5 for a careful discus- 
sion) to extrapolate to the mass-shell of the phys- 
ical states. We therefore use the zero-momentum 
chiral Ward identities. For the two-point Green 
functions of the relevant operators these are[|| 

2n/(Q Q}S a0 + Mac{<t>% Q) = 
2n/(Q 4)S a0 + AU^l cfe) + $ ah = 

(19) 

which imply 

M ac M bd {cpl 4) - -(M<5>) ab +(2n f ) 2 (Q Q)5 a0 6 bo (20) 

We also need the result for the general form of 
the topological susceptibility (see eq.(12)): 

-A 



x(0) = (Q Q) 



1 - (2n / )M(M$) 00 1 



(21) 



Although the pseudoscalar operators 0g and Q 
couple to the physical states ?y Q = (ir° , 77, 77'), it is 
more convenient to redefine linear combinations 
such that the resulting propagator matrix is di- 
agonal and properly normalised. So we define op- 
erators r\ a and G such that 

( (Q Q) (Q 4) \ f(GG) 











(22) 



{G G) = -A 

and we impose the normalisation 
-1 



-5 a0 



(26) 



(27) 



This implies that the constants f aa in (25), which 
are simply the decay constants, must satisfy the 
(Dashen) identity 

f aa ml f Tf3b =<M<&&»^*«» 



-(M$) q6 + (2n f ) 2 A6aoS bo 



(28) 



The last line follows from the Ward identities (20) 
and (21). In terms of these new operators, the 
anomaly equation (14) at zero momentum is: 



/ 



2nfG5 a0 



(29) 



The notation — * is to emphasise that the iden- 
tity is only true for insertions of the operators 
into zero-momentum Green functions and matrix 
elements. At non-zero momentum, other opera- 
tors appear on the r.h.s. (In particular, we can 
not take the on-shell matrix elements between 
the vacuum and the \rj a } and use (0\G\r] a ) = 
to conclude that the decay constants f aa can 
be identified as usual with the matrix elements 
(0|5 M J^ 5 \r] a ) . As we have repeatedly emphasised, 
the anomaly removes the familiar correspondence 
between the decay constants and the current ma- 
trix elements.) 

Now recall how conventional PCAC is applied 
to the calculation of 7r° — > 77. The pion decay 
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constant is denned as the coupling of the pion to 
the axial current 

(01^)=^ => (0\d»J* 5 \n) = Uml (30) 

and satisfies the usual Dashen formula The next 
step is to define a 'phcnomcnological pion field' ir 
by 



(31) 



To include electromagnetism, the full anomaly 
equation is extended as in (5) to include the 
F^F^ contribution. Using (29) we have 

^(77l^ 3 5 |0) 

= Uml <77k|0) + (771^^10) 

= f n ml(n n)( 71 \n) + a^^F^F^O) 



(32) 



where (tt tt) is the pion propagator — l/(fc 2 — m\\ 
At zero momentum, the l.h.s. vanishes because of 
the explicit factor and the absence of massless 
poles. We therefore find, 



3 " 
a cm~ 
TT 



(33) 



In the full theory including the flavour singlet 
sector and the gluonic anomaly we find a similar 
result. The 'phcnomcnological fields' are defined 
by (29) where the decay constants satisfy the gen- 
eralised Dashen formula (28) . Notice again that 
they are not simply related to the couplings to the 
axial current as in (31) for the flavour non-singlet. 
We therefore find: 

i^(77l^ 5 |0)=/ aa m2 /3 (77|^|0) 

+2n / <77|G|0>* a0 + < m ^( 17 \F^F^\0) 

+2n f (G G)( 11 \G)S aO +a- m ^(j 1 \F^F^\0) 



using the fact that the propagators are diagonal 
in the basis ?y Q , G. Using the explicit expressions 
(26), (27) for the Green functions, we find in this 
case: 



where the extra coupling gc-yy is defined through 
(34). This completes the derivation. It is evi- 
dently a straightforward generalisation of conven- 
tional PCAC with the necessary modification of 
the usual formulae to take account of the extra 
gluonic contribution to the axial anomaly in the 
flavour singlet channel, the key point being the 
identification of the operators 77" and G in (29). 

3. Phenomenology 

In this section, we discuss critically the way in 
which data on 77' (77) — > 77 decays (and closely 
related processes such as 7/(77) ~~ * ^7) are anal- 
ysed in the phenomenological literature^. The ex- 
isting analyses are based on formulae in which 
the impact of the colour C/a(1) anomaly has not 
been correctly taken into account. We therefore 
propose a comprehensive re-analysis of the data 
based on the formulae (8), (9) derived above. 

The two-photon decay widths are given by 



Y{t)'{t}) -> 77) = 



64tt 



I sy (77)77! 



(36) 



The current experimental data, quoted in the 
Particle Data Group tables [O are 



r(r7->77) = 0.510 ± 0.026 keV 



(37) 



arising principally from the 1988 Crystal Ball 
and 1990 ASP @ results on e+e~ -> e + e~r] 
(here, we follow the note in the 1994 PDG com- 
pilation and use only the two-photon 77 pro- 
duction data), and 



r(77' -> 77) = 4.28 ± 0.19 keV 



(38) 



dominated by the 1998 L3 data jl5j on the two- 
photon formation of 77' in e + e~ — > e + e~ ir + n~ 7 . 

For the purposes of comparing our theoretical 
results with the standard phenomenology, we fo- 
(34) cus on t ne very thorough and complete discussion 
in rcf.jL6|, updated in (Tt]]. These authors have 
studied not only the two-photon decays of 77' and 



f aa 9v a ii + 2rifA (7G77 $a0 



a a - 



(35) 



3 Here, and in section 4, we have only cited the few pa- 
pers which have been most influential in the preparation 
of this article. This is not intended to be a comprehen- 
sive review of the diverse and interesting literature on the 
phenomenology of r\ and rf physics. 
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r\ but also the related radiative vector-meson de- 
cays 77' (77) — > Vj, where V = (p, w, </>). They also 
analyse ip — > 77' (77)7 decays. 

There are two main points of difference between 
our theory and the analysis in these papers (and 
all other phenomenological treatments of which 
we are aware, including the notes on the 77' and 77 
decay constants and two-photon decay formulae 
in the Particle Data Group tables pi]] ). The first 
is that the term involving the coupling <?g77 m 
the decay formula (8) is neglected. The second is 
that the decay constants entering the decay for- 
mulae are assumed to be given (just as for the 
pion) by the matrix elements of the axial cur- 
rent. Further, in [|l6| it is assumed that the decay 
constant matrix (10) has just three independent 
components, although this is generalised to four 
in the second paper Q . 

To make this concrete, let us define quantities 

faa by 

<0|d% a 5 |7 7 ' 3 ) = rml p (39) 

The decay formulae used in |l6],|l7| are then sim- 
ply 



r a 9v> 



(40) 



in our notation. The most obvious problem with 
eq.(40) is that it is not consistent with the renor- 
malisation group. Since the singlet axial current 
J° 5 is not RG invariant because of the anomaly 
(see section 5) neither are the f 0a . This is recog- 
nised in |l^,|l^l but the associated running of the 
'decay constants' f 0a is assumed there to be neg- 
ligible. However, this is not the main problem, 
which is that the RG non-invariance of f 0a is a 
signal that the basic decay formula (40) is wrong 
- as we have seen above, it must be modified to 
include the coupling gc-y-y- 

The interesting question for phenomenology to 
address is whether 5G77 is in fact small, as sug- 
gested by the OZI-based discussion above. The 
success of the programme carried through in 
IG.Tt]] suggests that this may well be true. How- 
ever, the only real way to settle the question is to 
repeat the analysis using the correct formulae (8) 
and to fit the data with the full set of decay con- 
stants f aa (subject to the Dashen constraints (9)) 



and the extra coupling gc^-y- The actual value of 
gc-ry may also be of relevance in other related 
processes (see below). 

The second main problem stems from the mis- 
taken identification of the relevant decay con- 
stants with the axial current matrix elements 
(39). It is clear that the anomaly equation may 
be used to relate the (non RG invariant) matrix 
elements (0|Q|t7 q ) of the topological charge den- 
sity to the 'decay constants' f aa . Specifically, 
if we neglect m u ,md, we have in our notation 
(ciM0) 



{0\Q\rf 



l a0 



V3f 8a ml 



(41) 



This expresses the ratio of the coupling of Q = 
fjtrG^G^ to 77' and 77 in terms of the f°v'M. 
It has been pointed out in 18 that the matrix 
elements (0|<3|r7'(?7)) enter into the formulae for 
the radiative decays ip — > 77' (77)7. The current 
data fill is 

T(V> -> 77 7 )/r tota i = (0.86 ± 0.08) x 10- 3 (42) 
and 

r(V> -> r/7)/r t otai = (4.31 ± 0.30) x 10~ 3 (43) 
where 



total 



(87 ± 5) keV 



(44) 



Studying these decays therefore gives information 
on the 'glue content' of the 77' and 77. However, in 
refs.[16 Dj] the formula (41) is used in conjunction 
with the erroneous (40) to relate the (0\Q\ri' (ri)) 
to the decay constants appearing in 77' (77) — > 77. 
This is formally incorrect, and in any case it must 
surely be inconsistent to use the naive decay for- 
mula (40) , which is derived on the basis that both 
77 and 77' are true pseudo-Goldstone bosons, to es- 
timate their anomaly-induced glue content. In 
truth, as originally pointed out in [jj, the physi- 
cal, RG invariant, decay constants f aa appearing 
in the true decay formulae for the 77' and 77 are in 
principle quite independent of the quantities f aa 
related to the axial current matrix elements. 

A final comment on the theory concerns 77 — 77' 
mixing and the decay constants f aa . For some 
time, it was usual to describe 77' (77) — > 77 decays 



8 



in terms of two decay constants (/° and / 8 ) and 
a single 77 — 77' mixing angle defined in terms of 
the SU (3) eigenstates by 

\r] ) = cos?%o) + sin$|fy 8 ) 
\v) = - sin 1 770) + costf|?7 8 ) 

(45) 

More recently [^9|- p2| , ^7| etc., it has been realised 
that the f aa decay constant matrix (10) should 
instead have four independent components, which 
are conventionally parametrised in terms of the 
two / 0,s and two angles. In the version favoured 
by |L7| , these would be defined through 



Z 8 "' / 8 "J 



7°cos^ -/°sin^\ 
v / 8 sin^ / 8 cos^ 



Results for and § n derived using (40) are 
quoted in ref.|17|], where the difference between 
the two angles is claimed to be an 'energy depen- 
dence' of the 77 — rf mixing angle, with one be- 
ing evaluated on the 77' mass-shell and the other 
on the rj mass-shell. To us, however, it seems 
that this interpretation has no consequence. It is 
clearly correct to use four parameters to describe 
the decay constant matrix f aa and a parametri- 
sation in terms of two constants and two angles 
is as good as any other, but we see no reason 
to superpose the notion of an 'energy-dependent 
mixing angle' on this simple formalism. 

The analysis of rf{rf) — > 77 decays presented 
here can clearly be extended to study 77' (77) — » 
V7, where V is a flavour-singlet vector meson 
p,tj,(p, using the familiar assumptions of vector 
meson dominance. This allows us to extract the 
couplings g v '(ri)Vj from the 3-point Green func- 
tions (0\J°f J™ 1 J^|0) by relating them to the 
appropriate electromagnetic anomaly coefficients 
in exactly the same way as for the two-photon de- 
cays, then assuming pole-dominance for the vec- 
tor current. This is explained in detail in ref . |Tg| . 
However, exactly the same comments apply. The 
true decay formulae for P — > V7 involve extra 
gav-y couplings in the singlet sector (the full set of 
formulae analogous to those quoted in ref. |]l6| can 
be trivially derived using the method described in 
section 2). These arc omitted in ref. 16 . Again, 
the apparently reasonable agreement with data 



[fllfl found there suggests that these new couplings 
may be small, but this should be checked by re- 
peating the analysis including the gav-y and find- 
ing the best fit to data. It would also be interest- 
ing to make a detailed comparison with the quark 
model/OZI analysis of ref. f23[| . 

Another obvious extension is to the decays 
77' (77) — * 7r + 7r~7, related to the 4-point Green 
functions (0|j"jf J+ 5 J~ 5 J™|0). This Green 
function has an AAAV box anomaly in direct 
analogy to the AVV triangle anomaly consid- 
ered above. Clearly this can be analysed sim- 
ply using the Ua(X) PCAC theoretical meth- 
ods described here. On the experimental side, 
while 77 — > 7r + 7r _ 7 has been observed with decay 
widthlBUl 



r(»j -> 7r+7T~7) = 0.056 ± 0.005 keV 



(47) 



the most recent L3 results find no evidence for 
the non-resonant 77' — > 7r + 7r~7 decay in a channel 
dominated by the p: rf — > fr/, p — > 7r + 7r~. This 
is in contradiction with earlier claims J24|-^6| that 
the non-resonant process had been seen. 

To conclude this discussion of radiative 77' and 
77 decays, it will be clear that a re-analysis of the 
experimental data using the decay and Dashcn 
formulae (8), (9) should be performed to verify 
whether or not gcy-y is negligibly small or, if not, 
to measure it. We would then have a theoreti- 
cally reliable determination of all the parameters 
playing a role in this sector. This would be free 
of any additional prior theoretical input based on 
either OZI |||| or the large N c expansion in the 
chiral Lagrangian framework [ plp^j27| -|29|] . At 
the same time, it should be recognised that the 
'glue content' (0\Q\rj' (i])) of the 77' and 77 are not 
expressible in terms of the true, RG-invariant de- 
cay constants f aa and the relations (41) should 
not be used in the analysis of ip — > 77' (77)7 data. 

4. Related experiments 

We now mention briefly a few other experi- 
ments where the theoretical techniques described 
in this paper can be applied. 

The first is the behaviour of the polarised pho- 
ton structure function gj, considered as a func- 
tion of the 'target photon' momentum p. Stan- 



9 



dard analysis of the two-photon DIS process 
e + e _ — > e + e~X reduces the problem of find- 
ing the first moment of gj to a non-perturbative 
evaluation of the AVV three-point Green func- 
tion (0\J° 5 Jl m (p)J c p m (-p)\0). This is of course 
precisely the problem in evaluating 77' (77) — > 77, 
except that we now want to study a range of 
off-shell photon momenta. This has been car- 
ried out in rcfs. . A sum rule for gj has 
been presented which fixes the p — ► and asymp- 
totically large p behaviour of the first moment 
of gl in terms of the electromagnetic and colour 
Ua{1) anomaly coefficients. For intermediate val- 
ues of p r~j trip, the behaviour is shown to be sen- 
sitive to the non-perturbative realisation of chiral 
symmetry breaking. Detailed results are given in 
refs.gfl). 

The theoretical techniques developed here are 
equally applicable to the flavour-singlet contribu- 
tion to the first moment of the polarised proton 
structure function g^, the famous 'proton spin' 
problem. This is shown to be related to 
a Ua{1) Goldberger-Treiman relation, which is 
evaluated easily by the techniques of section 2 and 
which involves an analogous new coupling gcNN- 
The OZI and RG conditions in this case are how- 
ever rather different |?3|-|36f . 

More recently, it has been proposed f3?],|38| that 
the theoretical formalism of section 2 may also 
provide a way of understanding recent data from 
the CLAS collaboration 39 4(J on the photo- and 
electro-production of <j> mesons. The idea here is 
that one important contribution to the photopro- 
duction process 7 A — * <j)N may be modelled by 
'glue' exchange, which in turn could be related 
to the goNN coupling in the J7a(1) Goldberger- 
Treiman formula and the gc^-y coupling extracted 
from 7/(77) — > (fry as discussed above. This is an 
interesting proposal which certainly deserves to 
be pursued. 

Another closely related process is 7/ (77) photo- 
production. Both 7/j — * rjp plj ] and 773 — > r/'p 
|42|]43jl have been studied in low energy experi- 
ments close to threshold. A theoretical discus- 
sion of these reactions, based on an effective ac- 
tion incorporating the UaO-) PCAC formalism of 
section 2, may be found in ref.Q]. For a recent 
general review of rf physics along these lines, em- 



phasising the importance of the J7a(1) anomaly 
and gluonic degrees of freedom and incorporating 
an extensive survey of experiments, see also 45 
and references therein. 

We should also make a special mention here 
of the new WAS A 4tt detector @ at CELSIUS, 
which will provide an important facility for pre- 
cision 77 and 77' decays. This is of course the stim- 
ulus for the present Workshop. 

Finally, and somewhat remarkably, the process 
77' (77) — > 77 plays an important role in calcula- 
tions of the strong-interaction contributions to 
the anomalous magnetic moment g^ — 2 of the 
muon. This is of great current interest because of 
the discrepancy between the most recent exper- 
iments [ fl7| , f48| and theory (see e.g. Q for a re- 
cent compilation of sub-process contributions to 
g p — 2) and the suggestion that this could be a 
signal of non-standard model physics (e.g. p0|]). 
The important contribution for our purposes is 
the 0(a 3 ) correction arising from the Feynman 
diagram in which a QCD sub-diagram is linked 
to the muon line by three virtual photons, ef- 
fectively a hadronic contribution to light-by-light 
scattering. The dominant light pseudoscalar in- 
termediate state therefore involves back-to-back 
P — > 7*7* vertices, where of course the photons 
are off-shell. This is fully discussed in rcf. |Ti] 
where it is recognised that the conventional theo- 
retical analysis is inadequate and a proper treat- 
ment incorporating the gluonic anomaly along 
the lines discussed here is really required. In 
its absence, ref. |3l"| relies on using phenomeno- 
logical models and inputs from experiment. The 
result is that the total hadronic light-by-light 
scattering contribution to the muon g^ — 2 is 
-79.2(15.4) x 10" 11 , to be compared with the 
current experiment versus theory discrepancy in 



= !Cff„-2) of, 



,SM _ 



333(173) x 10" 



5. Theory 

In this section we derive the 77' (77) — * 77 de- 
cay formula using the method of operator Green 
functions and generalised 1PI vertices developed 
in refs. jl],|j etc. We present the discussion in three 
parts - first the anomalous chiral Ward identities 
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for the Green functions of the relevant composite 
operators and the 1PI vertices, then the deriva- 
tion of the decay formulae in terms of 1PI vertices 
(which are essentially the couplings gy (77,(3)77 )i 
and finally the renormalisation group equations. 
We then comment briefly on alternative methods 
involving chiral Lagrangians. 

5.1. Chiral Ward Identities 

The anomalous chiral Ward identities for QCD 
with massive quarks have been written down in 
the form used here in the review ||. We refer 
to this article for more complete derivations and 
restrict ourselves here to the most essential iden- 
tities. For the moment, we omit the electromag- 
netic contribution to the anomaly. 

The composite operators involved in the Green 
functions and 1PI vertices studied here are the 
currents and pseudoscalar operators J° 5 , Q, (j)^ 
and the corresponding scalar (f> a . We use the com- 
pact notation for the quark masses and conden- 
sates introduced at the start of section 2. 

The expressions for these operators given in 
section 1 define the bare operators. The renor- 
malised composite operators arc defined as fol- 
lows: 

tO _ 7 tO ja^Q _ tQ#0 

J fj,5 — ZjJ t ihB "V5 — J /i5B 

Q =Qb-^-zWj% b 

<i>i =z^ B r = z^b 

(48) 

where Z^ is the inverse of the mass renormalisa- 
tion, Zrj, — Z~ . The non-trivial renormalisation 
of J° 5 means that its matrix elements scale with 
an anomalous dimension 7 related to Z . This oc- 
curs because J° 5 is not a conserved current, due 
to the anomaly Q. Notice in particular the mixing 
of the operator Q with <9 M J° 5 under renormalisa- 
tion. As explained in |52j| , this leaves the combi- 
nation (9 M J" 5 — 2nfQ) appearing in the anomaly 
equation invariant under renormalisation. 

The Green functions for these operators are 
constructed by functional differentiation from the 
generating functional W{V® 5 ,9, , S a ], where 
V^ 5 , 9, S§ , S a are the sources for the composite 
operators J^ 5 ,Q,4>^,4> a respectively. For exam- 
ple, the Green function i(0\T Q(x) Q(y)\Q) is 



given by se ^)Ye( y ) > wnidl we abbreviate as Wee- 
This compact notation is perhaps unfamiliar, but 
is very convenient for manipulating complicated 
expressions involving the chiral Ward identities. 

In this functional formalism, the anomalous 
chiral Ward identity is 

d^Wv^ = 2n f 5 aQ We + M ac W SS 

- d adc S d W st + d adc StW S c 

(49) 

This can be compared with eq.(15). Notice the 
presence of the source variation terms on the r.h.s. 
This immediately give the identities for the 2- 
point Green functions: 

ikfjWy^yb^ - 2n f S a oWe V b 5 - M ac W S a V ^ = 
ikpWv.o ~ 2n f 8 aO Wg - M ac W s *e = 
ik^Wya^b - 2nf5 a0 W eS b - M ac W s * S b - $ ab = 

Combining these individual equations, we find the 
familiar identity 

k^kvWy*^ - M ac ^ cb = W S a DS t o (51) 

where Sfj is the source for the current divergence 
operator D a = 2n,fS a oQ + M ac §\. 

For our purposes, we really only need the zero- 
momentum chiral Ward identities. Clearly as- 
suming there is no massless boson coupling to the 
L0t(l) current, these are just 

2n f S a0 Wgg + M ac Wsie = 
2n f 5 a0 W es <, + M ac W S c Sb + $ a6 = 

(52) 

which implies the following identity for the topo- 
logical susceptibility 

(2n / ) 2 X (0) = M 0c W sss «M d0 + (M$) 00 ( 53) 

The 1PI vertices used below are defined as func- 
tional derivatives of a second generating func- 
tional (effective action) T, constructed from W 
by a partial Legendre transform with respect to 
the fields Q, and (j> a only (not the currents 
J° 5 ). The resulting vertices are '1PF w.r.t. the 
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propagators for these composite operators only. 
This separates off the particle poles in these prop- 
agators, and gives the closest identification of the 
field theoretic vertices with the physical couplings 
such as g n <* 77 ||. 

The basic anomalous chiral Ward identity for 
r follows immediately from that for W: 

dfjT V o. & = 2n f S a0 Q + M ac 4>l 



(54) 



and other identities follow simply by functional 
differentiation. In particular, for the two-point 
vertices, we find the following identities 



*aJV» B = 



iknTya.Q - 2n f S a + ®acF$°Q = 

ik^Tya^b + $ OC r c 6 - M ab = 

from which follows 



kf_i k v Yya yb 



cb 



(55) 
(56) 



At zero momentum, these reduce to 
®acT<t>cQ ~ 27^/(5 a o = 

^acT^b - M ab = 

which together imply 



(57) 



(58) 



The fact that the topological susceptibility is 
zero for vanishing quark mass can be seen imme- 
diately from eq.(53). One of the simplest ways to 
derive the precise form (12) or (21) is in fact to 
use an identity involving T. The two-point ver- 
tices are simply the inverse of the two-point Green 
functions (propagators), so in the pseudoscalar 
sector we have the matrix inversion formula: 

T QQ = -{W ee - Wes^WsssXbWs^y 1 
= -{Wee - W eSs M ac (MWs 5 s,M)^M db W S b e 



and using the identities (52) and (53) this implies 



r Q ^ = -x(i-(2n / ) 2 x(Mi>)oo 1 



(60) 



all at zero momentum. Inverting this relation 
gives the important result for the topological sus- 
ceptibility, 



X=-r Q Ul-{2n f fT- Q UM^) 



)oO 



(61) 



Substituting the explicit expression for (M$) 00 1 
(which is easily found from the definitions above), 



( 2n f) 2 ^ m q (qq} 



(62) 



we see that (61) reproduces the general form 
(12) where we can now identify the (mass- 
independent) non-perturbative coefficient as 



A ~ T QQ 



(63) 



We have already exploited these formulae in sec- 
tion 2. 

5.2. 7/(17) — > 77 from 1PI Vertices 

We are now ready to present the derivation 
of the decay formula (3) and generalised Dashcn 
formula (4). The technique relies on the iden- 
tification of the couplings g^yy with the zero- 
momentum limit of the appropriate 1PI vertex 
functions introduced above. 

The starting point is the Ward identity (54) 
extended to include the electromagnetic contri- 
bution to the anomaly for the axial current: 



d^T V a 5 = 2n f 8 aQ Q + a a em Q exa {A) + M ac (f> c 5 

— dacd^^tpi + dacdCjjfr^a 



(64) 



Q C m(A) is just shorthand notation for ^F^F^" , 
where F^ u is the field strength for the electro- 
magnetic field A^. (Since we are working only to 
leading order in a, it is not necessary to consider 
Qem as an independent composite operator with 
(59)non-trivial renormalisation.) 
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Differentiating twice w.r.t. the field A^, evalu- 
ating at the VEVs, and taking the Fourier trans- 
form, we find 



V%A x Ap 



r, a a c m<*ivfi 

" a em- e A<r Q /3Pl P 2 
TT 

— d a bc4> ^4, b A x AP 



(65) 



where p\ , pi are the momenta of the photons. No- 
tice that the mass term in (64) does not con- 
tribute explicitly to this formula. From its def- 
inition as 1PI w.r.t. the pseudoscalar fields, the 
vertex Tya^ A \ AP has no pole at k 2 — (even in 
the chiral limit) so the first term vanishes at zero 
momentum k, leaving simply 



'.A X AP 



k=0 



a 

i 

7T 



(66) 



(To simplify notation, it will be convenient 
from now on to define vertices T with the 
kinematical factors removed, e.g. T,pa A x AP = 

~r^a A x AP t\ aa pv\V2 ■ ) 

The first step in converting (66) to the de- 
cay formula (3) is to identify the physical states 
rf- . These appear as poles in the propagator ma- 
trix for the four pseudoscalar operators Q, ^ 
(a = 0,3,8). To isolate these poles, we diago- 
nalise the propagator matrix in this sector then 
normalise the three operators coupling to the 
physical states. 

We therefore define the operator 



G = Q-W dS «(W S5 s 5 )- b 1 4 



(67) 



so that by construction the propagators (G </>§) 
all vanish. (Notice that integrations over re- 
peated spacetime arguments are implied in this 
condensed notation.) Then define operators 



7f = C ab (t>\ 

such that the propagator matrix 

( V a = w sas „ = c aa w s?sb c Tb e 



(68) 



77 ^7) 
1 



/ k 2 -m 2 



V o 





-1 

k' 2 -m' 2 
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(69) 



where are the sources for the operators -q a . 

This change of variable affects the partial func- 
tional derivatives in T,f,a A \ AP in (66), which in- 
volves j^vt at fixed Q. In terms of the new vari- 
ables G, r] a we have 

_ 5ri a 8 SG 6 
q ~ 5<j)% Sr] a + S^5G 



6 S 
- - (Ws 5 s 5 )-t w s*ejQ 



Sij 



(70) 



The decay formula therefore becomes 



$abC Tba T vaA x Ap - <bab{Ws 5 S 5 ) a bW S b e T GA x A „ 



— "cm „ 
TT 



(71) 



The decay constants are identified as 
f aa = <S> ab C Tba 



(72) 



In terms of the propagators, we can write (from 
eq.(68)) 



f aa (Ws n sXef TI3b = *ac(Ws B S B )Zi*db 
and so at zero momentum 

t a rnl P f TPb = *ac{Ws B s B )£*db 



(73) 



(74) 



as quoted in (28). 

The remaining steps in finding (3) and (4) are 
an exercise in manipulating the zero-momentum 
Ward identities. First note that combining the 
two identities in (52) gives 



M ac W S c Sd M db 



-(M$) a6 +(2n / ) 2 X (0)<5 aO (5 6o (75) 



whose a, b — component is just (53). Note that 
(M$) a b is symmetric. Also define loo = S a0 6bo- 
Then we can write 



= ($M) ac (MWs B s B M)^M de W SS e 
= -2n/(M*) oc (-(M$) + (2n/) 2 X (0) 

= 2n /X (0)(l - (2n / ) 2 X (0)(M$) oo 1 ) _1 <5 aO 

= - 2n f T QQ S a0 



I-.k.) X X(0) 

cO 



(76) 



13 



where in the final step we have used the identifi- 
cation (61). Similarly, 

= (<S>M) ac (MW S5 s 5 M)^(M<S>) db 

= (Af$) ac (-(Af$) + (2n / ) 2 x(0)l O o)" 1 (M$) db 

V / cd 

= -{M<^) ab + {2n f ) 2 Y Q 1 Q 5 a0 5 b0 



(77) 



f aa F v c, A >. A p + 2nfT Q Q T GA x A „ — a" m — (78) 



This establishes the required results. Substi- 
tuting (76) and (77) into (71) and (74) we find 
the decay formula 

a 

i 

7T 

where the decay constants satisfy 

f aa <pf Tt3b = -(M^) ab +(2n f ) 2 T Q 1 Q S a0 S b0 (79) 

The final step is to identify the 1PI vertices 
with the couplings defined in section 1, viz. 



- n a A x AP — 9v a 71 



(80) 



and similarly for ga-yy- It is at this point that 
the central dynamical assumption is made. In 
fact, eqs.(78) and (79) are exact identities, fol- 
lowing simply from the definitions and the zero- 
momentum chiral Ward identities. To make con- 
tact with the radiative decays of the physical par- 
ticles, we must assume in particular that the 1PI 
vertex evaluated at k = accurately approxi- 
mates the physical coupling^, which is defined 
on mass- shell. This requires that T n a A \ AP has 
only a weak momentum dependence in the range 
< k 2 < m 2 a . This is reasonable, since it is 
defined to be a pole-free, amputated, RG invari- 
ant dynamical quantity. However, as in standard 
PCAC, the assumption is expected to be excellent 
for the 7r but progressively worse as the mass of 
the pseudo-Goldstone bosons increases. The hope 
here, in common with all attempts to include the 
T]' in the framework of PCAC (including chiral 
Lagrangians with 1/N C effects included |H|,p2 27 



|29[ , is that the approximation remains sufficiently 
good at the mass of the r( . 

4 The assumption that the 1PI vertices as defined here can 
be identified at all with the decay couplings of the physical 
particles rests on pole dominance - the assumption that 
the dominant particle poles in the pseudoscalar propagator 
matrix are indeed those of the r} a (see eq.(69)). 



5.3. Renormalisation Group 

It is important to determine the renormalisa- 
tion group behaviour of all the quantities appear- 
ing in these formulae. Recall, for example, that 
the RG behaviour was a crucial factor in the con- 
jecture that gG-y-y may be small in the leading 
OZI approximation. In general, the RG equa- 
tions play a key role in understanding the physics 
of the U A (1) channel. We therefore include here 
a brief and rather novel discussion of the RGEs 
for the relevant Green functions and 1PI vertices 
in the functional formalism. The essential results 
were first given in ref . Q , but are generalised here 
to include SU(3) breaking and r) — rj' mixing. The 
results are a straightforward extension of refs.jl],[3| 
but were not explicitly written down in 

The fundamental RGE for the generating func- 
tional W in pure QCD follows immediately from 
the definitions (48) of the renormalised composite 
operators. It is: 



VW = i V ^-2n- f 



Wy 



s a w s * + 



(81) 



where T> 



The notation + . . . (which is suppressed in the fol- 
lowing equations) refers to the additional terms 
of 0(k 2 ) and 0(k A ) which are required to pro- 
duce the contact term contributions to the RGEs 
for n-point Green functions of composite opera- 
tors. (This notation is omitted in the following 
equations, but it should be remembered that it 
is implicit.) These terms are discussed fully in 
refs.|53 lj, but will be omitted here for simplicity. 
They vanish at zero momentum so do not directly 
affect the derivation of the decay formulae, but 
do have important implications for the validity 
of PCAC extrapolations from zero-momentum to 
on-shell quantities. 

The RGEs for Green functions are found simply 
by differentiating eq.(81) w.r.t. the sources. Sim- 
plifying the results using the chiral Ward identi- 
ties (50) , we find a complete set of RGEs for the 
2-point functions. These are: 



VWy 



(7<5a0 



l5 0b )W V a V b 
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VW v %e = hS aa + i)W V \9 + l7^—M ob W V a s „ 

M5 /»• 2.11 r 

VWya S b = (7(5 a O + J^Wya gb 

11.R R 1 1. R R 



1 



PWftsj = (7 + 7*)W" fls » + T^jCMteWsjsj + $ ob ) 



The pattern of cancellations which ensures the 
consistency of these equations with the chiral 
Ward identities is quite intricate, but may readily 
be checked. 

At zero momentum, we can immediately use 
the second of eqs.(52) to write the above RGE 
for W$e as 



VWee = 2lW 9d 

fc->0 



- 2 ^77T-^( M oaW S a Sb M bo 



(M$) c 



= 



(83) 



using (53). This shows that the zero- momentum 
topological susceptibility is a RG invariant, 



V X (0) = 
and thus 
VA = 



(84) 



(85) 



where A is the non-perturbative parameter in 
eq.(12), which enters into the final decay and 
Dashen formulae. 

Next, we need the RGE for the generating func- 
tional of the 1PI vertices. This follows immedi- 
ately from its definition and the RGE (81) for W: 



f)T = "f(v® 5 — — !— 

"70 (Vs 



b a T, 



(86) 



where V = I U-i-+B-§-— r y m ^air— ' 

Y'd^^dg lrn^ q q dnlq y y Q ^^ 

The RGEs for the 1PI vertices are found by dif- 
ferentiation, and using the Ward identities (55) to 



simplify the results, we find for the pseudoscalar 
sector: 

VT QQ = -2 7 r QQ + 2 7 -^- hocTQQT^c 
2n f L 

VT Q4,l = -(7 + 7^)r Q0 6 



+7 



^ r ^$ = -2 70 r 0? 
l 



+7 



2n f 



Here, 2? = 2? + 7 0(0 Q ) 3 | r . As explained in rcf.|§, 
this is identical to the RG operator V defined 
above (acting on W) when the sources are set to 
zero and the fields to their VEVs. 

These RGEs play two roles in the discussion. 
First, they are used as consistency checks on the 
various formulae we derive. Second, according to 
our conjecture, they provide the clue to identi- 
fying quantities which are likely to show viola- 
tions of the OZI rule and those for which we may 
reasonably expect the OZI limit to be a good 
approximation. This is because we can identify 
quantities which will be particularly sensitive to 
the Ua(1) anomaly as those which have RGEs 
involving the anomalous dimension 7 . 

We now derive the RGEs for the Green func- 
tions and 1PI vertices involved in the various ex- 
pressions related to the 77' (77) — * 77 amplitude. 
To do this, we first need to include the electro- 
magnetic fields and their anomalous dimensions. 
As already noted in section 5.2, the anomalous 
dimension for the composite operator Q C m(A) en- 
tering the anomaly equation is of 0(a 2 ) so can be 
neglected at the order at which we are working. 
We denote the anomalous dimension correspond- 
ing to the usual electromagnetic field renormali- 
sation by 7^4. 

The RGEs for the 1PI vertices Tq A x Ap an d 
T^b A \ AP are easily found by differentiating 
eq.(86) and simplifying using the Ward identities. 
We find, 



(V + 2x4) T QA x AP = 
-lT QA x AP +7- — T QQ a° n - e\ paS 3ViV2 

ZUf IT 
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and, 



(88) 



(V + 2 lA ) T^ MAP 
1 



l<tFd>lA*AP +17, «em- ^\pa(iVlV2 

° ZTlf 5 7T 



1 



These are very similar to the corresponding equa- 
tions in ref. ffl , with the obvious inclusion of the 
mass term and SU (3) breaking in the VEVs and 
flavour mixing. 

These expressions simplify remarkably at k 2 = 
0. Using the zero-momentum chiral Ward identi- 
ties (57) for Tq^o and T^c^b, together with (66), 
we find 



(V + 2j a )T qa x Ap 
{V + 27 j4 )r b j4 A j4 p 



fc=0 ' ' ^ 6 

From the latter, we immediately have 



= 

fc=0 

-70I> A a Ap (O) (90) 



(V + 2 lA )$> ab T^ A x AP 







fc=0 



(91) 



verifying the RG consistency of the basic identity 
(66). 

It only remains to rewrite these results in terms 
of the 1PI vertices for r\ a and G. First, recall the 
identification of the decay constants: 

Tba 



where 



(92) 



(93) 



The l.h.s. is the propagator matrix for the phys- 
ical rj a and is therefore RG invariant. From the 
RGE for W S a S b, we then find 



VC aa = -7 C Qa 

and it follows immediately that 

Vf aa = 



(94) 



(95) 



This confirms that the true decay constants f aa 
are RG invariant. Contrast this with the current 



matrix element definition (0|J° 5 |?7 Q ) = ik^f aa , 
for which 



(96) 



we imme- 



Vf aa = lS a0 f aa 

Now, since T GA x AP = i QA x AP , 
v" ' <t>t 
diately deduce from eq.(90) above that 

(V + 2 lA ) 9Gll = (97) 

Finally, from eq.(70), we have 
(89) 

^ab^^M AP — f aa ^ri a A x AP 

-^ab{Ws,S 5 Y a lW sig T GA x AP 

(98) 

Combining (90) with the RG identities (82), we 
then find after further use of the chiral Ward iden- 
tity (52) that indeed 



(V + 2 7A )T vaMAP -0 

k=0 

at zero momentum, i.e. 
(V + 2 lA )g, ril = 



(99) 



(100) 



as promised. In fact, if we had included the con- 
tact terms in the RG equations throughout, as in 
ref. jl) , we would have found at this point that the 
coupling 5r)°77(fc 2 ) is actually not RG invariant 
for all k. However, it was found in |]J] by keep- 
ing careful track of the contact terms that it is 
also RG invariant on- shell. This is an important 
point - it is a necessary condition for the dynami- 
cal assumption that the on-shell couplings may be 
well approximated by their zero-momentum val- 
ues, essential to the PCAC method, to be valid. 

This completes our survey of the RG properties 
of the radiative rf{rj) — > 77 decay formulae con- 
firming that, as stated in section 1, all the quanti- 
ties appearing in the formulae are RG invariant. 
In particular, this confirms the identification of 
f aa as the true, physical decay constant. 

5.4. Chiral Lagrangians, OZI and 1/N C 

An alternative to the approach presented here 
is the popular method of chiral Lagrangians, so 
we include a few comments on their relation. Chi- 
ral Lagrangians are models of low-energy QCD in 
which the basic fields are chosen to parametrise 



1G 



the coset manifold G/H (for a chiral symmetry 
breaking pattern G — > i?) and thus lie in one-to- 
one correspondence with the Goldstone bosons. 
The dynamics, which is determined by the isom- 
etry group of this coset manifold, is therefore 
arranged from the outset to satisfy the (zero- 
momentum) chiral Ward identities. The great 
advantage of chiral Lagrangians is that they pro- 
vide a systematic way of going beyond leading 
order in a low-momentum expansion, higher or- 
der terms being developed by the loop expansion 
in this non-renormalisable QFT fl54| . 

It is, however, important not to forget that chi- 
ral Lagrangians are simply models of QCD. They 
implement the chiral Ward identities in a particu- 
larly elegant, geometric way but they still implic- 
itly assume the same dynamical approximations 
of pole dominance (in selecting the most relevant 
low-energy states) and smoothness of momentum 
extrapolations that are explicit in the actual QCD 
treatment in terms of operator Green functions. 

This is especially important when chiral La- 
grangians are extended |2l| , |2^ , |27| -|29[l to the non- 
linear U(3) x U(3)/U(3) models incorporating the 
77', which is of course not a Goldstone boson be- 
cause of the anomaly. The dynamics of these 
models is therefore not entirely constrained by the 
geometry of the coset space but must be imple- 
mented in part by hand if they are to be accurate 
representations of true QCD. The most promis- 
ing systematic approach is to use the 1 /N c expan- 
sion, since at leading order in 1 /N c the 77' becomes 
a true Goldstone boson (because the anomaly is 
sub- leading) . However, as we have emphasised, 
the leading 1/N C (or the OZI []) approximation, 

5 In the text, we have preferred to refer to the OZI approx- 
imation, rather than l/N c . The OZI limit is precisely de- 
fined H as the truncation of full QCD in which non-planar 
and quark-loop diagrams are retained, but diagrams in 
which the external currents are attached to distinct quark 
loops, so that there are purely gluonic intermediate states, 
are omitted. (This last fact makes the connection with 
the familiar phenomenological form of the OZI, or Zweig, 
rule.) This is a more accurate approximation to full QCD 
than either the leading 1/N C limit (N c — * 00 at fixed rif), 
the quenched approximation (n f — > at fixed N c ) , or the 
leading topological expansion (N c — > 00 at fixed nj/N c ). 
In the OZI or leading 1/N C limits, the Ua(X) anomaly is 
absent, there is an extra Goldstone boson, and there is no 
meson-glueball mixing. 



while a good approximation for some quantities, 
is completely invalid for others. In our presen- 
tation, we have pursued the consequences of the 
anomalous chiral Ward identities as far as possi- 
ble without making extra dynamical assumptions, 
introducing these only at the end to make contact 
with the physically observed couplings and decay 
constants. In particular, we have used the renor- 
malisation group as a guide to which quantities 
we expect to have a smooth 1/N C perturbation 
expansion and which violate 1/N C or OZI signifi- 
cantly at leading order. 

It would therefore be of considerable interest 
to make a detailed comparison of the 1 /N c chiral 
Lagrangian predictions pl| , p2p^ , |29| ] with those 
made here (and also for the closely related analy- 
sis of the J7a(1) Goldberger-Treiman relation p3| - 
|36[ ] and its link with the 'proton spin' structure 
function g\). Since the fundamental anomalous 
symmetry and dynamical assumptions should be 
the same, it would be interesting to see how these 
are realised in these two, in principle equivalent, 
approaches. 

6. Epilogue 

In this paper, we have reviewed the theory and 
phenomenology of the radiative 77' (77) — » 77 de- 
cays, together with closely related processes such 
as rf(rj) -> V7, 7/(77) ~> tttTj i> v'ivh, 
etc. The theory of these decays is indeed a 
tale of two anomalies: first, the electromagnetic 
?7a(1) anomaly, which was spectacularly success- 
ful historically in explaining the otherwise mys- 
terious 7r — ► 77 decay; second, the gluonic £^a(1) 
anomaly, which makes the physics of the flavour- 
singlet 0~ channel in QCD so subtle and interest- 
ing. 

Indeed, it is the role of the gluonic J7a(1) 
anomaly that makes the 77' and its decays worth 
studying. It is therefore disappointing that this 
new physics is so often obscured by phenomeno- 
logical analyses which try to fit the data into the 
straightjacket of decay formulae written down in 
naive analogy with tt — > 77, without taking the 
implications of the gluonic anomaly fully into ac- 
count. The purpose of this paper is to urge a 
fresh phenomenological look at 77' physics, treat- 
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ing both the electromagnetic and colour Ua(1) 
anomalies in a complete and theoretically self- 
consistent manner. 

The goal, going beyond mere confirmation of 
the well-established physics of pseudo-Goldstone 
bosons and their interactions, is to gain new phe- 
nomenological insight into the rich and fascinat- 
ing subject of gluon topology in QCD. 
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